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ABSTRACT 

We investigate the nature of the substantial population of high-redshift galaxies with J s -K s > 2.3 colours 
recently discovered as part of our Faint InfraRed Extragalactic Survey (FIRES). This colour cut efficiently 
isolates galaxies at z > 2 with red rest-frame optical colors ("Distant Red Galaxies," or "DRGs"). We select 
J s -K s > 2.3 objects in both FIRES fields, the w 2.5' x 2.5' Hubble Deep Field South (HDF-S) and the «5'x5' 
field around the MS 1054-03 cluster at z = 0.83; the surface densities at ZQ.Vega < 21 mag are 1.6 ±0.6 and 
1.0 ±0.2 arcmin" 2 , respectively. We here discuss a sub-sample of 34 DRGs at 2 < z < 3.5: 11 at X^vega < 
22.5 mag in HDF-S and 23 at ^T s .vega < 21.7 mag in the MS 1054-03 field. This sample enables for the first 
time a robust assessment of the population properties of DRGs. We analyze the A = 0.3-2.2 /im spectral energy 
distributions (SEDs) constructed from our very deep near-infrared (NIR) and optical imaging collected at the 
ESO Very Large Telescope and from the Hubble Space Telescope. We develop diagnostics involving the I$i4- 
7 S , J s -H, and H-K s colours to argue that the red NIR colours of our DRG sample cannot be attributed solely to 
interstellar dust extinction and require for many the presence of an evolved stellar population with a prominent 
Balmer/4000 A break. In the rest- frame, the optical colours of DRGs fall within the envelope of normal nearby 
galaxies and the ultraviolet colours suggest a wide range in star formation activity and/or extinction. This 
is in stark contrast with the much bluer and more uniform SEDs of Lyman-break galaxies (LBGs). From 
evolutionary synthesis modeling assuming constant star formation (CSF), we derive for the DRGs old ages, 
large extinctions, and high stellar masses, mass-to-light ratios, and star formation rates. For solar metallicity, a 
Salpeter IMF between 0.1 and 100 M Q , and the Calzetti et al. extinction law, the median values for the HDF-S 
(MS 1054-03 field) sample are 1.7 (2.0) Gyr, A v = 2.7 (2.4) mag, M* = 0.8 (1.6) x 10 11 M , MjL v>i , = 1.2 
(2.3) MqLy'q, and SFR= 120(170) M Q yr _1 . Models assuming exponentially declining SFRs withe-folding 
timescales in the range 10 Myr- 1 Gyr generally imply younger ages, lower extinction, and lower SFRs, but 
similar stellar masses within a factor of two. Compared to LBGs at similar redshifts and rest-frame V-band 
luminosities, DRGs are older, more massive, and more obscured for any given star formation history. For the 
entire sample of ,ff s -band selected galaxies in the FIRES fields at 2 < z < 3.5 and to the above magnitude limits, 
we find that the derived ages, extinctions, and stellar masses increase with redder J s -K s colours. Although the 
rest-frame optical colours of the DRGs are similar to those of local normal galaxies, the derived properties are 
quite different; detailed studies of this new z > 2 population may significantly enhance our understanding of 
how massive galaxies assembled their stellar mass. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — infrared: galaxies 



1. INTRODUCTION 

Understanding how galaxies formed and evolved is a cen- 
tral challenge of modern astronomy. In the past decade, spec- 
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tacular progress in instrumentation and observing techniques 
has enabled major advances as the high-redshift universe has 
been opened up for research. Our current view has been 
largely influenced by the discovery of an abundant population 
of actively star-forming galaxies at redshi fts z ~ 3 -4 selected 
by th e efficient Lyman-break technique JSteidel & Hamilton! 
119931 Steidel et al. 1996a b). These Lyman-break galaxies 
(LBGs) are among the best-studied classes of high-redshift 
objects: large samples have been sp ectroscopically con- 
firmed at z ~ 3-4 (e.g. ISteidel et al. IH999t) and extensive 
investigations have focussed notably on their stellar popu- 
lations, star formation histories, che mical abundances, kine- 
matics, and clustering properties (e.g.lGiaval isco et alTI|1998t 
i Pettini et al. 1 1200 It iPapovich. Dick inson. & Ferguson 1200 It 
Shanl ev et al. Tl200lll2003UFrb et al. II2003I) . LBGs dominate 
the ultraviolet (UV) luminosity density at high redshift and 
their number density is comparable to that of L* galaxies lo- 
cally, making them a major constituent of the early universe. 
In the context of the widely favoured hierarchical scenarios 
of galaxy formation, they are thought to be the progenitors 
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of present- day massive galaxies in groups and clusters (e.g. 
Baugh etal. 1199a. 

Yet, by construction the Lyman-break technique relies on 
a strong Lyman discontinuity in the rest-frame far-UV and 
is necessarily biased towards relatively unobscured galaxies 
with intense recent or on-going star formation activity. The 
typical stellar ages derived for LBGs are of a few 10 8 yr 
with moderate extinction Ay < 1 mag (Parjo vichetal]l200ll 
Shanl evet al. l l2001). Older and more quiescent systems at 
Z ~ 3 that have formed the bulk of their stars at z > 4 or, al- 
ternatively, more obscured galaxies may have esc aped detec- 
tion in optical surveys (e.g. Ferguson. Dickinson. & Panovich 
12001 . With the advent of 8-10 m-class telescopes and 
the development of sensitive near-infrared (NIR) instruments 
equipped with large-format detectors, it has become possible 
to extend very deep surveys to longer wavelengths and ac- 
cess the rest-frame optical emission of sources at z ~ 1-4. 
Compared to the rest-frame UV, the rest-frame optical spec- 
tral energy distribution (SED) of galaxies is less affected by 
the light of young massive stars and by dust extinction, and 
better traces longer-lived stars that dominate the stellar mass. 

Colour criteria involving NIR bandpasses can now be 
applied to identify new populations at z > 1. The 
long-standing debate over the formation of massive early- 
type galaxies has driven considerable interest in the dis- 
covery of, and subsequent systematic searches for, red 
objects such as the "extremely red objects" (EROs) 
generally defi ned by R — K > 5 or I - K > 4 mag 
(e.g. iFlston. Rieke. & Riekel [19881 IhTi & Ridgwavl 11991 
Thom pson et al. 1119991 |Yan et al. 1 120001: IScodeggio & Silval 
2000[ iDaddi et al. ll2000t lMcCarth v et al. 1120011 ISmith et al. I 
2002, among others). In our own Faint InfraRed Extra- 
galactic Survey (FIRES; iFranx et al.ll2000l) . based on very 
deep NIR imaging of the Hubble Deep Field South (HDF- 
S) and of the field around the z = 0.83 cluster MS 1054-03, 
we have identified significant numbers of fairly bright (K s ps 
19 — 22.5 mag) candidate z > 2 galaxies selected from their 
Js~Ks > 2.3 colours. Analysis of the HDF-S sample sug- 
gests that this population makes a comparable contribution 
to the stellar mass dens ity at z ~ 3 as LBGs ( Rudnic ket al. I 
2003; Franx et al. 2003) and hence may be a substantial com- 
ponent in terms of stellar mass. Remarkably, there are far 
fewe r such bright red objects per u nit area in the HDF-North 
(e.g. lDickinson et al. 120001 see also lLabbe et al. 120031) . Cos- 
mic variance may however largely explain these differences 
since both Hubble Deep Fields are small and the J*-Kp > 2. 3 
population may be strongly clustered dDaddi et al. I 2003). 
Candidate high-redshift galaxies with unusually red J s - K s 
colours have been reported by other authors as well (e.g . 
Scodeggio & Silval2000tlHall et al. l200lllTotani et al. I200U 
Saracco et al. 2001, 2003). However, the focus has been on 
those objects with the most extreme colours, which are rarer 
and mostly at the faintest magnitudes. 

An immediate question is what causes the red colours of 
these objects. Very red NIR colours can be produced by an 
evolved population at z > 2 due to the stellar photospheric 
Balmer/4000A break redshifted into the J band and beyond 
or by high levels of extinction in galaxies possibly at lower 
redshifts. The presence of a highly obscure d active galac- 
tic nucleus (AGN) is another possibility (e.g.. lNorman et ai~l 
2002; Koekem oeret al. 120041) . While it appears that a combi- 
nation of stellar population aging and extinction effect s is re- 
quired to explain the observed colours (e.g. lDickinson et al. I 



2000tlTotani et al. l2001llHall et al. 1200 fl IFranx et al. 12001 
Sarac coet al. 1 12003). the results so far remain inconclusive 
due to the scarcity of systematic studies of large samples and 
to insufficient observational constraints. 

In this paper, we take advantage of our FIRES data set to 
address the above issue from the A = 0.3-2.2 /xm broad-band 
SEDs of 34 objects with J s -K s > 2.3 selected from both fields 
surveyed. The addition of the MS 1054-03 objects triples the 
original HDF-S sample. We focus on the ensemble pro perties 
of the J s -K s red galaxies. Ivan Dokkum et al. I J2003LI2004T) 
present results of our follow-up optical and NIR spectroscopy 
of a bright sub-sample in the MS 1054-03 field, which are 
particularly relevant to this work. 

We briefly describe the data in §|2] We discuss the selection 
criteria applied to construct the Js-K s selected samples in §[3] 
In § we analyze the properties of the SEDs to investigate 
the stellar populations and dust obscuration. We model the 
SEDs using evolutionary synthesis in §|5]to constrain quan- 
titatively the stellar ages and masses, the extinction, and the 
star formation rates. In §[6] we interpret our results and com- 
pare the derived properties to those of LBGs and of Zf-band- 
selected objects at similar redshifts in the FIRES fields. We 
summarize the paper and main implications in §0 For conve- 
nience, we adop t the term "distant red galaxies," or "DRGs," 
introduced by van Dokku m et al. 1 J2001 to designate candi- 
date z > 2 galaxies with J s -K s > 2.3 colours. All magni- 
tudes are expressed in the Vega-based photometric syste m ex- 
cept w hen "AB" indicates reference to the AB system (Oke 
119711) . Throughout, we assume a A-dominated cosmology 
with fl m = 0.3, il A = 0.7, and H = 70 h 7Q kms" 1 Mpc" 1 . 

2. DATA 

The data were obtained as part of FIRES iFranx et al. I 
2000), a public NIR survey of the HDF-S and MS 1054- 
03 fields carried out at the European Southern Observatory 
(ESO) Very Large Telescope (VLT). The observations, data 
reduc tion, and sou rce catalogues are described in detail by 
Labbe et al .1(12003) for HDF-S and N.M. Forster Schreiber et 
al. (2004, in preparation) for the MS 1054-03 field 9 . 

Briefly, the NIR observations were taken in the 7 S , H, and K s 
bands with the VLT Infrared Spectrograph And Array Camera 
(ISAAC: lMo"orwood et al. 11991 . For HDF-S, a total of 103 h 
integration time was spent in a single « 2.5' x 2.5' pointing 
covering the Hubble Space Telescope (HST) WFPC2 main 
field. We complemented the NIR data with the publicly avail- 
able deep optical WFPC2 imaging in the t/300, B450, V606, and 
/gu bands ICasertan oet al. 120001) . For the MS 1054-03 field, 
77 h of ISAAC integration time was obtained in a w 5' x 5' 
mosaic of four pointings. At optical wavelengths, we used 
existin g HST WFPC2 m osaics in the V$o6 an d ^8 14 bands 
l lvan Dokkum et aT. 1 1200(1) and collected additional imaging 
in the Bessel U, B, and V bands with the VLT FORS1 in- 
strument. With HDF-S, FIRES achieves the currently deep- 
est ground-based NIR imaging and the deepest /f-band map 
to date, even from space. With the MS 1054-03 field, the 
area surveyed is nearly five times wider, down to ss 0.7 mag 
brighter magnitudes. 

The sources were detected i n the band using v ersion 
2.2.2 of the SExtractor software ( Bertin & Arnouts 1996). For 
consistent photometry across all bands, the fluxes were mea- 
sured on the maps convolved to a common spatial resolution 

9 The reduced images, photometric catalogues, and photomet- 
ric redshifts are available online through the FIRES homepage at 

|http : / /www . strw . leidenuniv . nl/~f ires| 
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of full-width at half maximum FWHM = 0.48" for HDF-S 
and 0.69" for the MS 1054- 03 field, matching to the map 
of poorest seeing in each data set (H and U band, respec- 
tively). From the curve-of-growth analysis of the average stel- 
lar profile in each band, the fractional enclosed flux agrees to 
within 3% at diameters 0.7" < d < 2.0" for HDF-S and 2% at 
1.0" < d < 2.0" for the MS 1054-03 field, where the diam- 
eters range from the smallest (1.5 x FWHM of the matched 
PSFs) to the largest apertures relevant for the colour mea- 
surements. Colours and SEDs used in this work are based 
on measurements in custom isophotal apertures defined from 
the detection map, hereafter referred to as "colour" fluxes or 
magnitudes. Total magnitudes in the K s band were computed 
in apertures based on auto-scaling apertures (Rron 1980) for 
isolated sources and adapted isophotal apertures for blended 
sources. The photometric uncertainties were derived empir- 
ically from simulations on the maps. The total 5er limiting 
magnitudes for point sources are = 23.8 mag for HDF- 
S and 23.1 mag for the MS 1054-03 field. We will refer to 
the objects by their number in catalogue versions v3.0e for 
HDF-S and v3.1b for the MS 1054-03 field, prefixed with 
"HDFS-" and "MS-." 

Wherever appropriate, we corrected for the Galactic extinc- 
tion along the HD F-S and MS 1054-03 lines of sight ( A v = 
0.09 and 0.1 1 mag: ISchlegel. Finkbeiner. & Davisll998l) . We 
also accounted for the gravitational lensing by the foreground 
MS 1054-03 cluster based on the lensi ng model constructed 
by iHoekstra. Franx. & Kuiikenl (|2000) from weak lensing 
analysis. The average magnification factor for all Zf s -band se- 
lected objects considered in this paper is 1.18 (or 0.18 mag). 
We emphasize that the colours and SED shapes are unaffected 
since lensing is achromatic. 

Photometric redshifts z p h were determined by applying an 
algorithm involving linear combinations of redshifted spec- 
tral templates of galaxies o f various types, as described by 
iRudnick et al. I J2001L l2?X)3h . Monte-Carlo simulations were 
used to estimate the z p h errors accounting for uncertainties 
in the fluxes as well as template mismatch, and reflecting 
not only the confidence interval around the best solution but 
also the presence of secondary solutions. Comparison with 
available spectroscopic redshifts z sp imply an accuracy of 
5z= <|z S p-Zph |/(1 +z sp )_) =0.074 for both fields. Forthez > 2 
regime of interest in this work (see § !3.2t . Sz = 0.051 based 
on the 20 /f s -band selected sources with z sp determinations. 

3. NEAR-INFRARED SELECTED RED Z > 2 SAMPLES 

3.1. The J s — K s Colour Criterion 

To identify red candidate high-redshift galaxies, or DRGs, 
we applied the cr iterion J s - K s > 2.3 mag discussed by 
Fran x et al. I §003). This criterion was specifically designed 
to target the stellar photospheric Balmer and 4000 A breaks 
at redshifts above 2, which are the most prominent features 
in the rest-frame optical continuum SED of galaxies. The 
Balmer discontinuity at 3650 A is strongest in A-type stars 
and, more generally, is indicative of stellar ages ~ 10 8 - 
10 9 yr. At older ages, the 4000 A break characteristic of 
cooler stars with types later than about GO, and strongest in 
giants and supergiants, becomes dominant. It is due to the 
sudden onset of metallic and molecular opacity bluewards of 
4000 A produced notably by Ca II H+K, Fe I, Mg I, and CN 
lines. 

Based on synthe t ic spe ctral templates computed with the 
iBruzual & Charloi d2003l) models (see § |5j, unextincted 



single-age stellar populations older than w 250 Myr or pas- 
sively evolving galaxies formed at Zf > 5 satisfy the criterion 
at 2 < z < 4, where the Balmer/4000A break is redshifted 
into the NIR regime. Galaxies forming stars at a constant 
rate have a J s -K s colour redder than 2.3 mag only if they 
are > 1 Gyr old and moderately obscured at z > 2, or younger 
and highly obscured with Ay > 2 mag at z > 1 . As argued by 
van Dokkum et al. 1 J200l l2004). our follow-up spectroscopy 
is consistent with the J s -K s > 2.3 criterion selecting z > 2 
galaxies with a high efficiency of ~ 80%. 

The J s - K s colour cut at 2.3 mag is slightly blue r than 
used in other recent studies (in particular T otani et al. 11200 1^ 
Sarac coet al. 12 003). Our objective is not to focus on the rare 
extreme objects but to allow inclusion of systems with rest- 
frame colours characteristic of normal nearby galaxies, which 
may be missed by rest-frame UV-selected surveys. Observed 
colours of J s - K s > 2.3 mag correspond to rest-frame colours 
of B-V > 0.5 or U-B > 0.1 mag at z « 2.5. These rest- 
frame colours encompass those of the majority of present-day 
luminous gal axies, from elliptical to inte rmediate-type spiral 
galaxies (e.g. Buta et al . 1 1994t Hansen et al. 120001) . 

As discussed bv !Labb6eHil*T d2003l) . contamination by 
foreground cool stars is expected to be negligible since most 
known M, L, and T dwarfs have 7 S - K s colours bluer than 
2.3 mag; others such as extreme carbon stars or Mira variables 
that could be redder are rare and unlikely to be present in pen- 
cil beam survey s at high ga lactic latitudes such as our FIRES 
fields (see also lHaU et aUl2001tlSaracco et al. 1120031) . From 
NIR spectroscopy of six DRGs in the MS 1054-03 field and 
the Ch andra Deep Field South (CDF-S), Ivan Dokkum et al. I 
(2004) argued that strong emission lines in the K s band are 
unlikely to account for the observed red colours of most of 
the population. The main contaminants are more likely very 
dusty (Ay > 2 mag) galaxies at 1 < z < 2, probably at the 
~ 20% level (Ein Dokkum et al. 120031) . 

Another potential concern is that samples selected with the 
Js~K s > 2.3 criterion suffer from contamination by AGN. 
However, the majority of the bright DRG s we identified in 
the HDF-S and MS 1054-03 fields (§E3 are well resolved 
in our A^-band imaging data, implying that their NIR broad- 
band fluxes are not strongly influenced by continuum emis- 
sion from an AGN. Emission lines from type II AGN could 
also influence the NI R colours . The rest-frame optical spec- 
tra of van Dokkum et al. ( 2004) show that this effect is minor 
for the confirmed type II AGN in our samples, and we are 
pursuing a NIR spectroscopy program to confirm this result 
10 . For the moment, we conclude it is unlikely that emission 
from AGN affects the NIR fluxes and colours in a substantial 
way. 



3.2. Sample Selection 

For the purpose of SED analysis and modeling, our goal 
was to construct a sample of J s -K s selected candidate z > 2 
galaxies with robust NIR colours, reliable SEDs, and span- 
ning a redshift range that would allow consistent comparisons 
with the well-studied LBGs. We required that the J s -K s > 
2.3 mag objects have a signal-to-noise ratio S /N > 15 on their 
,?if s -band colour flux and a minimum of 30% of the total expo- 

10 The fairly high fraction of ~ 40% of type II AGN in the spectroscopic 
sample i van Dokkum et al. 2003 , 2004 ) is likely to result from the spectro- 
scopic selection bias towards galaxies that are bright in the observer's optical 
and with bright emission lines. 
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sure time (or relative weight w > 0.3) in all bands 11 . This 
yielded initial samples of 14 objects at Kl ot < 22.5 mag in 
HDF-S and 31 at < 21.7 mag in the MS 1054-03 field. 
The effective sky-plane areas are 4.48 and 23.86 arcmin 2 , re- 
spectively. 

The photometric redshifts of the selected sources lie in the 
range 1.88-4.26. Spectroscopic redshi fts are available for 
four of them in the MS 1054-03 field Jvan Dokkum et al. I 
I2003l [2004: S. Wuyts et al. 2004, in preparation). Three have 
z sp between 2.42 and 2.61, and their z p h agrees with z sp within 
1 a. The fourth has z p h = 1 -94^ ^ and turned out to be an inter- 
loper at z sp = 1 . 1 89. We note that the large uncertainty towards 
low redshifts for this source reflects the presence of a long tail 
of Zph solutions with comparable likelihood. Figure [2 com- 
pares the Zph's with z sp 's and shows the redshift distributions 
of the samples, where we used z sp whenever available instead 
of Zph- For the present study, we focused on the J s -K s selected 
sources with redshift between 2 and 3.5. We chose this inter- 
val because we are primarily interested in galaxies at z > 2 
and it matches well the range for LBGs selected from space- 
based HST WFPC2 data. Ground-based selected LBGs typ- 
ically lie at somewhat higher redshift with less overlap with 
the range covered by our DRG samples. This is mainly be- 
cause the selection criteria rely on a different set of filters 
with , notably, a redder "U" ban dpass (for a discussion, see, 
e.g..lGiavalisco & Dickinsonl200ll). 

The samples finally adopted, with J s -K s > 2.3 and 2 < z < 
3.5, consist of 11 objects at Kg* < 22.5 mag in HDF-S and 
23 at fC° l < 21.7 mag in the MS 1054-03 field. The red- 
shift distribution for the HDF-S (MS 1054-03 field) sample 
has a mean (z) = 2.60 (2.44), median (z) med = 2.50 (2.42), 
and dispersion er(z) = 0.42 (0.30). These values are com- 
puted using the z sp whenever available, and are essentially 
unchanged when considering only the z p h's (the largest dif- 
ference is (zph)med = 2.30 for the MS 1054-03 field). There 
are four objects just below the z = 2 limit whose la z p h error 
would allow a z > 2 and, conversely, five objects just above 
whose la z p h error would allow z < 2; we verified that these 
small sample variations would not alter the conclusions of our 
SED analysis and modeling. Table^gives relevant character- 
istics of the adopted DRG samples. 

We note that not all the MS 1054-03 field objects for which 
spectroscopy is presented by Ivan Dokkum et al. H20031I2004I) 
are included in our sample. The initial selection for the spec- 
troscopic follow-up was based on a preliminary version of 
the reduced data and photometry. In the final photometric 
catalogue (v3.1b), the two galaxies with spectral signatures 
indicating an AGN are slightly bluer than the J s - K s cut- 
off (MS -1035 and MS- 1356, numbered 1195 and 1458 in 
Ivan Dokkum et al. l l2003). The bright = 19.77 mag and 
red J,-K s = 2.69 ±0.09 m ag object MS - 140 (number 1 84 in 
Ivan Dokkum et alH 120031) satisfies all criteria except for the 
minimum integration time in the NIR bands, with relative 
weights w(J s ) = 0.15, w(H) = 0.21, and w(K s ) = 0.18. How- 
ever, it has a spectroscopic redshift (included in Figure [Qj) 
and, at z sp = 2.705, lies in the range considered for our anal- 

1 ' The weights in each band for the HDF-S are normalized to the respec- 
tive total exposure time over the entire field covered by the maps. For the 
MS 1054-03 mosaic, the NIR ISAAC fields have somewhat different total in- 
tegration times and the highest overall weights are reached in the small over- 
lapping regions. To ensure that the selection area covers fairly the mosaic, 
we considered the weights for the NIR data normalized to the total exposure 
time in each of the ISAAC fields individually. For sources in the overlapping 
regions, we required that w > 0.3 in at least one of the field. 



ysis. MS - 140 constitutes a valuable addition that we will 
discuss where relevant but keep out of the formal sample. 

Figure|2]shows the raw cumulative number counts per unit 
area as a function of total K s magnitude of the DRGs, for both 
cases of z > 2 and of no redshift restriction. To assess the vari- 
ations of surface density with limiting magnitude and between 
the two fields, the counts are extended to the 5er K^° x limits 
for point sources of 23.8 mag in HDF-S and 23.1 mag in the 
MS 1054-03 field. We did not account here for the lensing 
by the MS 1054-03 cluster. The corrections would be small, 
with an average dimming by « 0.2 mag and a mean reduction 
by » 20% from image- to source-plane area for the redshift 
range covered by the sources. These effects are opposite so 
that to first order, the surface density to a given apparent mag- 
nitude is insensitive to lensing (the exact effect depends on the 
slope of the luminosity function, which is unknown). 

Per unit area, the MS 1054- 03 field is richer in sources 
brighter than K^° x = 20 mag while HDF-S has higher surface 
densities for fainter limits down to 23 mag (where the com- 
pleteness levels still are > 90% for both fields). These differ- 
ences are not due to the lensing magnification by MS 1054- 
03. For z > 2, there are five DRGs at Kf l < 20 mag in the 
MS 1054-03 field and none in HDF-S; at fainter limits, HDF- 
S has higher surface densities by a factor of two on average. 
The redshift cut has little impact: without it, there are seven 
DRGs at Kl ot < 20 mag in the MS 1054- 03 field and still 
none in HDF-S, and the surface densities at fainter limits are 
1.7 times higher on average in HDF-S. The two fields cover 
small areas and cosmic variance could easily account for a 
large part of the differences in the surface densities, the more 
so in vi ew of the pos sible strong clustering of the DRG pop- 
ulation (Daddi et al. 2003). This might also explain the de- 
ficiency in HDF-North, with similar survey area as HDF-S, 
where we find no source with J s -K s > 2.3 colours at magni- 
tudes brighter t han A~ s = 21 mag using g round-based photom- 
etry (from Fernandez-Soto et al. 1999) 12 . 

Figure [5] plots the raw cumulative counts as a function 
of J s -K s cutoff for all K s -band selected sources at Kl ot < 
22.5 mag in HDF-S and < 21.7 mag in the MS 1054 -03 
field (the magnitude limits ensuring S/N > 15 on the /f s -band 
colour fluxes of the DRGs). For 2 < z < 3.5, the DRG sam- 
ples are not strongly sensitive to the exact colour cut: varying 
the criterion between 2.2 and 2.4 mag implies 1 I+3 objects 
for HDF-S and 23^ for the MS 1054-03 field, or 10% -30% 
variations. Imposing no redshift limits, the same colour cut 
variations result in 14^ objects for HDF-S and 3 1^| 5 for the 
MS 1054-03 field; the fraction of galaxies at z < 2 increases 
rapidly for J s -K s colours bluer than 2.3 mag. The median 
la uncertainty on the J s -K s colours among the adopted DRG 
samples is 0.13 mag and all those with uncertainties larger 
than 0.2 mag (the maximum is 0.43 mag) have colours redder 
than 3.0 mag. 

3.3. Comparison LBG Sample 

We also constructed a sample of LBGs taken in HDF-S. 
We used t he original HDF-S HST WFPC2 photometric cata- 
logue from Casertano et al. (2000) and selected the LBGs us- 
ing the criteria introduced bylGiavalisco & Dickinson ( 2001 ). 

12 We consider here specifically colours involving the ground-based J band 
because the corresponding F110W bandpass of NICMOS onboard HST is 
significantly bluer and about twice as broad as the J or J s bandpasses (the 
colour term is negligible between the latter two). Consequently, Fl 10W pho- 
tometry is not useful to select galaxies red in their / s — K s colours. 
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Of the 139 selected LBGs at V606.AB < 27 mag, 91 are cross- 
identified in our FIRES X" s -band selected catalogue and have 
w > 0.3 in our full set of maps. All LBGs in this sub-sample 
are brighter than Kl° l < 24.0 mag. To make up the comparison 
LBG sample, we then selected those satisfying the same A" s ! ot 
limit as for the HDF-S DRGs and with redshift 2 < z < 3.5. 
The requirement on the /T s -band magnitudes ensures that the 
selected LBGs have reliable NIR photometry and rest-frame 
V-band luminosities similar to those of the selected DRGs. 
The adopted sample consists of 33 LBGs at A" s tot < 22.5 mag. 
Their redshift distribution has (z) = 2.68, (z) me d = 2.76, and 
u(z) = 0.42. The mean and median redshifts are close to those 
of the 2 < z < 3.5 DRGs, especially from HDF-S. The aver- 
age (and median) J s -K s colour is 1.68 mag with dispersion 
a = 0.28 mag, and only two of these LBGs are redder than 
J s -K s = 2.3 mag. 

For this LBG sample, we will hereafter use the correspond- 
ing t/3oofi450 VatihwJsHKs photometry from the FIRES 
HDF-S catalogue (v3.0e). Although selected at different 
wavelengths, the use of SEDs from the same photometric cat- 
alogue ensures uniformity in the data sets for both DRGs and 
LBGs and more directly comparable results in the analysis 
and modeling. The requirement that our LBG sample be also 
detected in the K s band makes our selection slightly differ- 
ent than the classical purely optical selection. We did not 
construct a similar LBG sample for the MS 1054-03 field 
because of its ground-based UBV data, which would select 
LBGs in a somewhat higher redshift range than spanned by 
the DRG samples. 

4. PHOTOMETRIC PROPERTIES OF THE SAMPLES 

The large numbers of DRGs that we identified from / s - 
K s > 2.3 mag in the FIRES survey suggest they may represent 
a sign ificant component of the high-redshift gala xy popula- 
tion jFranx et al. II2003L Ivan Dokkum et al. l l2003). It is thus 
important to understand their nature, in particular their stellar 
populations and dust content. Before undertaking the detailed 
theoretical modeling, we here examine the broad-band prop- 
erties of the sample to assess the relative importance of aging 
and extinction in a qualitative but less model-dependent man- 
ner. 

4.1. Observed Optical to Near-infrared SEDs 

The observed SEDs of DRGs in both FIRES fields exhibit a 
wide variety in their shapes, especially in the optical regime. 
Figure0]presents the SEDs of selected representative objects 
in each of the HDF-S and MS 1054-03 field. The galaxies are 
sorted according to increasing Igi4-K s colour, and the SEDs 
are normalized to the flux density in the Ks band (the over- 
plotted SED fits will be discussed in § [5ji. At one extreme, 
some of the DRGs appear overall very red with little if any 
flux in the optical bands. At the other, some have blue ob- 
served optical colours indicative of on-going star formation 
activity. Most exhibit a rather abrupt transition redwards of 
the / s band reminiscent of a spectral break in the continuum. 
A few show a rather smooth and gradual increase in flux with 
wavelength suggestive of a featureless continuum reddened 
by dust obscuration. 

4.2. Rest-frame SEDs Properties 

The rest-frame optical colours of the DRGs lie within 
the range covered by present-day normal galaxies, in strik- 
ing contrast with the much bluer LBGs. Figure [5] shows 
the rest-frame SEDs of our samples of DRGs and LBGs at 



2 < z < 3.5. All SEDs are dereddened for Galactic extinc- 
tion. We accounted for the average intergalactic attenuation 
due to line blanketing by Lya ((Da)) and higher-order Lyman 
lines ((Db)) following the prescriptions of Madau ( 1995). The 
SEDs are compared with the empirical UV-optical template 
spectra of nearby galaxies of types E, Sbc, Scd, and Im from 
IColeman. Wu. & Weedmanl d!980ft . The spectra shown here 
were extended beyond their origin al A = 1400 — 10000 A cov- 
erage by Bolzonella, Miralles, & Pello (2000) using the evo- 
lutionary synthesis code of Bruzual & Chariot ( 1993, GIS- 
SEL98 version) with appropriate model assumptions. The 
SEDs and spectral templates are normalized to a flux density 
fx of unity at A rest = 2500 A using linear interpolation between 
the adjacent rest-frame data points. 

At A rest > 2500 A, the SEDs of all DRGs fall within the en- 
velope of normal galaxies. None of the DRGs has rest-frame 
optical colours as red as the earliest-type nearby ellipticals, 
perhaps not surprisingly so because the universe is only a few 
times 10 9 yr old at 2 < z < 3.5. At shorter wavelengths, the 
SEDs span a wide range in rest-frame UV slopes encompass- 
ing those of normal galaxies and extending to much redder 
slopes, suggesting a range in level of star formation activ- 
ity and in interstellar extinction. In contrast, the LBG SEDs 
exhibit less scatter especially in the UV and are systemat- 
ically bluer than lat e-type spirals, as first demonstrated by 
Papov ich et alJ ( 120011) . These authors showed that the SEDs 
of their HDF-N LBGs fit well in the envelope d efined by the 
local starburst templates of lKinnev et al. I ( 1199 6) and are con- 
sistent with the presence of a young and modestly obscured 
stellar population. 

Figure|6]illustrates the effect s of extinc t ion, w here we used 
the Im-type templat e of IColem an et al. | (11980) and applied 
the extinction law of Calzetti et al. ( 2000) assuming a simple 
uniform foreground screen of obscuring dust. The template 
SEDs inherently include some intrinsic reddening so that our 
simulations represent additional amounts of extinction. Inter- 
estingly enough, the Im template combined with added ex- 
tinction of Ay between and 3 mag brackets the range of 
rest-frame SED shapes of our DRG sample. Large extinction 
offers a natural interpretation for the many objects with very 
red rest-frame UV emission and consequently could also ac- 
count in part for the red rest-frame optical colours. Figures|5] 
and [6] do not allow us to assess the relative importance of 
evolved stellar populations and extinction effects in produc- 
ing the observed SEDs of DRGs. However, they demonstrate 
clearly that DRGs and LBGs are distinct in their rest-frame 
UV-optical SEDs, presumably because of very different dom- 
inant stellar populations and/or extinction levels. 

In Figure [5] the SEDs of the three spectroscopically- 
c onfir med DRGs among our sample along with MS - 140 (see 
§ I3.2> are plotted with different symbols. The comparison in- 
dicates that the galaxies in the spectroscopic sample do not 
deviate strongly in their rest-frame SEDs from the full pho- 
tometric sample. On the other hand, it also shows that the 
spectroscopic sample does not probe the entire range of prop- 
erties seen for the ensemble of DRGs. 

4.3. Optical Break versus Dust Extinction 

The wavelength sampling and resolution of our data are too 
coarse to delineate accurately a Balmer/4000 A break. Even 
when strongest, the break represents a moderate flux jump by 
a factor of ks 2, a subtle feature compared to the Lyman break 
at 912 A. For the brightest DRG of the spectroscopic sample, 
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the presence of a break is confirmed by the detection of a sig- 
nificant drop in continuum level in its NI R spectrum col lected 
with NIRSPEC at the Keck Telescope ( Ivan Dokkum et al~l 
2004). Here, we examine the observed NIR properties of our 
ensemble of DRGs using diagnostic colour-colour diagrams, 
which allow us to establish that a break is required to explain 
the colours of an important fraction of them. 

Specifically, we test whether the red observed NIR colours 
of DRGs (1) can be reconciled with a very young (< 10 8 yr) 
stellar population that is highly obscured by dust, or if (2) 
the presence of a more evolved (> 10 8 yr) stellar population 
whose SED has a prominent Balmer/4000 A break is neces- 
sary. In exploring various colour combinations, we relied on 
the sharpness of the Balmer/4000 A break versus the smoother 
wavelength dependence of dust attenuation. Extinction makes 
all colours redder while a prominent break produces a larger 
contrast between colours because of the steep drop in flux at 
the break. As a stellar population ages, the continuum slope 
intrinsically reddens because cooler stars dominate the opti- 
cal light but the break remains comparatively sharp. For red- 
shifted populations, these differences in colour contrast will 
be reflected in the variations with z in colour-colour space. 

We found lhaiH-K s versus J$-H and J s -H versus 1^4- J s 
provide the most sensitive diagnostic diagrams. The use of 
these adjacent bandpasses traces as closely as possible the 
Balmer/4000 A break in the range z ~ 1-3.5 while minimiz- 
ing reddening due to dust and to cool dominant stars. Fig- 
ure compares in these diagrams the colour distributions of 
the DRGs with the redshift evolution at 2 < z < 4 of se- 
lected model SEDs: 250 Myr and 1 Gyr old single-age stellar 
populations (SSPs), a passively evolving population formed 
at Zf = 5, and populations with constant star formation rate 
(CSF) and ages of 10 Myr, 100 Myr, and 1 Gyr. We com- 
puted the synthetic spect ra a nd colours using the codes of 
iBruzual & Charlotl (E0Q3b and lBolzonella et al. I (12000) as for 
the SED modeling and with the same standard assumptions 
(see §|5}. We also plot tracks for the 100 Myr old CSF model 
with Ay in the range 1-6 mag, applying the extinction law of 
iCalzetti et al. I (|2000) and a uniform foreground screen geom- 
etry. 

The model tracks in the top panels of FigureQshow the ex- 
pected trend that the older the age, the larger the amplitude of 
the colour variations with z as the stronger break moves red- 
wards across the successive bandpasses. The effect is maxi- 
mized for SSPs or passively evolving populations, where the 
break is not diluted by the light from younger stars with shal- 
lower or no absorption features and with bluer continuum. 
CSF models are consistently bluer in all colours and occupy 
narrower regions than SSPs of the same age. Extinction ef- 
fects shift the tracks along paths that run diagonally across the 
diagrams but preserve their amplitude at fixed age. As shown 
schematically in the bottom panels, the locus of CSF models 
with varying extinction gets wider with older ages, because of 
the stronger break. The effects of aging and extinction can be 
distinguished in that the loci of dusty model tracks broaden 
with increasing age in the direction perpendicular to the ex- 
tinction path, the key feature we here exploit. 

The immediate result from these diagrams is that the 
colours of the ensemble of DRGs cannot be reconciled with 
very young and highly obscured populations. Dusty CSF 
models need the broadening due to significant aging to bet- 
ter match the observed spread perpendicular to the extinction 
vector. To quantify the fraction of DRGs not overlapping 



with the dusty models, Figure [8] presents histograms where 
the colour distributions have been rebinned as a function of 
the deviation in / s — H and 7gi4 —Js around the line best repre- 
senting the dusty 100 Myr old CSF tracks (using robust linear 
regression). The maximum extent of the tracks about this line 
is indicated in the plots. The central loci and maximum widths 
for the 10 Myr and 1 Gyr old dusty CSF models are plotted 
as well. About two-thirds of the DRG sample has colours 
which are inconsistent with those of a 10 Myr old obscured 
star-forming population. The fraction drops to 30% -55% for 
an age of 100 Myr, and 5% -25% for 1 Gyr (depending on the 
colour considered). 

Interestingly, the H-K s versus J s -H diagram shows that if 
DRGs are actively star-forming, then the J s -K s > 2.3 colour 
cut will miss even those with a strong optical break unless 
they are highly obscured. In other words, the criterion pref- 
erentially selects not only for the most evolved star-forming 
systems but for the most obscured ones as well. The need 
for high extinction is alleviated for scenarios where the star 
formation rate is declining with time (as we discuss further in 
§|5j, but the J s -K s colour cut still selects for evolved galaxies. 
Figures Q and [8] demonstrate that the red NIR colours of the 
ensemble of DRGs cannot be attributed to extinction effects 
alone. For a significant fraction of the objects, the colours re- 
quire the presence of a strong Balmer/4000 A break produced 
by stars with ages > 250 Myr. 

5 . MODELING OF THE SPECTRAL ENERGY DISTRIBUTIONS 

In this section, we describe our model ingredients and fit- 
ting methodology. We present the results for our set of 
"standard" parameters with fixed initial mass function (IMF), 
metallicity, extinction law, and two choices of star forma- 
tion history (SFH). We also test the sensitivity of the results 
to changes in the SFH, metallicity, and extinction law. We 
modeled the broad-band SEDs of our DRG sample following 
the approach adopt ed in recent similar studies of high red- 
shift galaxies (e.g.lPapoyich et all200lUShaplev et al.ll2001b 
van Dokkum et al. 2004). We fitted redshifted synthetic spec- 
tra of stellar populations to the SEDs to constrain simulta- 
neously the age and extinction, keeping z fixed for each ob- 
ject. We derived the stellar masses (M + ) and star forma- 
tion rates (SFRs) by scaling those of the best-fit input model 
so as to reproduce the observed fluxes. We determined the 
stellar mass-to-light (M/L) ratios using the attenuated rest- 
frame V-band luminosities (Ly.*) computed as described by 
Rudnick et al. ( 2003, their Appendix C). Although the extrin- 
sic (i.e. dust-attenuated) luminosities from our best-fit mod- 
els agree to better than 10% o n average, we preferre d the less 
model-dependent estimates o flRudnick et al. I ( 120031) . 

5.1. Model Ingredients 

We generated the synthetic spectral templates with the 
most recent version of the evolutionary synthes is code de- 
velop ed by G. Bruzual and S. Chariot ( B ruzual & Charloll 
2003). This "BC03" version features improved stellar evo- 
lutionary tracks with observationally motivated prescriptions 
for the thermally-pulsing asymptotic giant branch phase. It 
also incorporates updated and newly available stellar libraries 
for two sets of spectral resolution and coverage. We used the 
lower resolution set covering A = 9 A to 160 /im relying on 
the BaSeL 3.1 library. We selected t he "Padova 1994" ey olu- 
tionary tracks, which are preferred by Bruzual & Chariot over 
the more recent "Padova 2000" tracks because the latter may 
be less reliable and predict a hotter red giant branch leading 
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to worse agreement with observed galaxy colours. We used 
the solar metallicity set of tracks. The only constraints avail- 
able for the chemical abundances of DRGs come from the 
[N II]A6584 A/Ha line ratio measured in two of the galaxies 
observed spectroscopicall y, which suggests h igh metallicities 
of Z ~ 1 - 1.5 Z Ivan Dokkum et al. 1120041 These DRGs 
appea r more metal- rich than the five LBGs at z ~ 3 studied by 
iPettini et al.lfeOOU 0. 1 -0.5 Z Q ) and similar to the seven UV- 
selected star-f orming "BX/MD" objects at z ~ 2 for which 
Shapl evet al. I J2004) inferred solar, and possibly super-solar, 
metallicities. In all cases, however, the determinations rely on 
limited samples and suffer from large uncertainties. 

The computed grid of spectral templates covers ages be- 
tween 10 5 and 2 x 10 ll) yr. The SFH most appropriate for 
our objects is unknown and broad-band p hoto metry alone 
is not effective at constraining it (see § 15.41 also, e.g., 
iPaoovich etaI]l200U [Shaplev et al. 1 120011) . We considered 
two star formation rate parametrizations: constant in time 
with R(t) = R(t = Q)=Rq and exponentially declining in time 
as R(t ) = Rq exp(-f jr) with timescale r = 300 Myr (hereafter 
"CSF" and "r3ooMyr" models). These exact choices are some- 
what arbitrary. The CSF model is intended to represent the 
case of significant on-going or recent star formation activ- 
ity. Evidence for this scenario in part of the DRG population 
includes the detection of Ha and/or [O III] AA4959,5007 A 
lines in the non-AGN so urces of the spectroscopic sample 
( Ivan Dokkum et al~l l2004) and the blue rest-frame UV SED 
of several others among the full sample (see Figures@]and[5}- 
The T3ooMyr model allows for the possibility of quiescent sys- 
tems that underwent a period of enhanced star formation in 
their past. 

We adopted a Salpeter IMF (dN / dm oc m~ 2,35 ) between 0. 1 
and 100 Mq. The IMF is unconstrained for our objects, and 
a steep rise down to the lower mass cutoff likely is unrealis- 
tic in view of the turnover below 1 M^ inf erred for the lo- 
cal IMF (e.g.. lKroupall200 It lChabrierll2 003 ). The most crit- 
ical impact is on the derived masses and mass-to-light ratios 
(M/L), which depend strongly on the shape and cutoff of the 
low-mass IMF. For instance, for a fixed rest-frame V-band lu- 
minosity, the derived masses an d M/L ra t ios wo uld be about 
a facto r of two lower with th elKrou pal d200 1 |) orlCh abrier 
120031) IMFs (e.g., Figur e 4 of Br uzual & Charlotll2003t see 
also IPapovich et al 1200 11 f or a discussion of the effects of dif- 
ferent IMFs in their SED modeling of LBGs). 

We have used the mass of stars still alive M t instead of 
the total mass of stars formed M tot . M tot is simply the inte- 
gral of the SFR and corresponds to the total gas mass con- 
sumed since the onset of star formation. is computed 
by subtracting from M tot the mass returned to the ISM by 
evolved stars via stellar w inds and supernova explosions (see 
iBruzual & Charlol2003l for details). We preferred M+ in our 
analysis because it represents the mass of the stars produc- 
ing the light. For the typical ages 1-2 Gyr of interest here 
and with the adopted IMF, the cumulative mass-loss effects 
are modest and comparable for the CSF and T3ooMyr models, 
reaching 20% -25% of M tot . Even in the limit of an SSP 
where the effects are largest, they amount to a similar frac- 
tion of 23% -26%. 

Both interstellar extinction by dust within the objects and 
attenuation due to intergalactic H opacity were applied to the 
BC03 templates while performing the fits. We explored a 
range of extinction with < Ay < 3 mag. Higher Ay values 
are physically possible but unlikely to be actually measured 



from rest-frame UV/optical data because the corresponding 
optical depths t\ =A\/ 1.086 largely exceed unity. The ex- 
tinction law and the geometry for the obscuring dust and emit- 
ting so urces are other unc ertain parameters. We adopted the 
law of iCalzetti et al. I d200 0) derived empirically from obser- 
vations of local UV-bright starburst galaxies under the formal- 
ism of a foreground screen of obscuring dust. As emphasized 
by these authors, this extinction law and model geometry may 
not be valid in dust-rich systems where the optical depths are 
large and the sources and dust are more or less mixed spa- 
tially. The average Lyman line absorption by the intergalactic 
m edium (IGM) w as accounted for based on the prescriptions 
of lMadauH 19951) . Lyman continuum absorption was approx- 
imated by setting the flux of the templates equal to zero at 
A res t<912A. 

5.2. Fitting Procedure 

We performed the fits using the publicly available HYPERZ 
photo metric redshift code, version 1.1 jBolzonella et al. I 
2000). We fixed the redshift to the z p h determined indepen- 
dently (see §|3 or to the z sp when available. The code was al- 
lowed to fit simultaneously the spectral template (age) and the 
extinction, for each of the SFH adopted. Interstellar extinction 
is first applied to the synthetic spectra, which are then red- 
shifted and further attenuated for IGM absorption. The fitting 
algorithm relies on least-squares minimization, with the chi- 
squared calculated as \ 2 = Z)J(/obs -bf te m P (z)) / a] 2 . The 
sum is performed over the n photometric bands. The template 
fluxes /temp(z) are computed from the reddened, redshifted, 
IGM-absorbed template spectra using the total effective trans- 
mission functions in each bandpass. The scaling factor b is 
determined from the full SED constructed from the colour 
fluxes, corrected by the same amount to the total aperture 
based on the Zf-band measurements; all SEDs have also been 
corrected for Galactic extinction and, for the MS 1054-03 
field sample, for the lensing magnification prior to the fitting 

(§□• 

We treated any flux with / D b s < 2a as a non-detection, with 
upper limit equal to the la uncertainty. The assigned la un- 
certainties and upper limits correspond to those of the flux 
measurements. We assumed a minimum error of 0.05 mag on 
the photometry to avoid the fits being driven by a few data 
points with very small errors. This also accounts for uncer- 
tainties in the absolute flux calibration. Our results are not 
strongly influenced by this assumption. If we do not impose 
a minimum photometric error, the median of the best-fit ex- 
tinction values remains the same for the two fields and two 
standard SFHs considered. The median M*'s vary by < 1%, 
and the median SFRs by < 15%. The median of the best-fit 
ages decreases from 1 .7 to 1 .0 Gyr for the HDF-S DRGs and 
CSF, and increases from 1.0 to 1.4 Gyr for the MS 1054-03 
field DRGs and the T300Myr model, not affecting our findings 
that DRGs harbour old stellar populations. 

HYPERZ resamples the age grid from the BC03 models 
from 221 to 51 fixed steps. We required that the template 
ages considered for the fitting do not exceed that of the uni- 
verse at the redshift of each source. We emphasize that the 
best-fit age, which we will denote f s f, is strictly speaking the 
time elapsed since the onset of star formation. An alterna- 
tive, perhaps more meaningful definition is the SFR-weighted 
mean age (t)sm. corresponding more closely to the age of 
the stars contributing the bulk of the stellar mass. For CSF 
models, (f)spR = 0.5 f s f. For exponentially decaying models, 
(Osfr = [?sf-T + rexp(-f s f/r)]/[l-exp(-f s f/r)], which be- 
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comes (f)sFR ~ fsf - 7" when f s f 3> r. Moreover, f s f should not 
be over-interpreted as the age of the entire galaxy. It is im- 
portant to remember that all the properties derived from the 
SED modeling apply only to the stellar population that domi- 
nates the fluxes in the observed wavelength range. This does 
not preclude, e.g., the presence of a pre-existing stellar pop- 
ulation t hat has faded out too much (a scenario explored for 
LBGs by Panov ich et alJ2 001 ) or of a younger population that 
is too obscured to be detected by our observations. 

We estimated the confidence intervals for individual ob- 
jects from Monte-Carlo simulations. For each source, we 
performed 500 synthetic realizations of the data by varying 
the fluxes by an amount randomly drawn from the (Gaussian) 
distribution of the measurement uncertainties (including the 
minimum error of 0.05 mag where appropriate). For sources 
with only a z p h determination, we also varied the redshift by 
drawing randomly from the associated probability distribu- 
tion P(z). We then repeated the fitting procedure for each re- 
alization. We determined the 68% confidence intervals from 
the distributions of best-fit values obtained with the simulated 
data. The behaviour of the confidence intervals for the DRG 
samples studied here are very similar to tho se of the spectro- 
scopic sam ple oflvan Dok kum e{ al~l (120041) and of the LBG s 
modeled bv lPapovich et all (1200 W and Shanl ev et aTl (120011) , 
i.e. often asymmetric around the best-fit values, reflecting 
strong degeneracies notably between f s f and Ay, and imply- 
ing thatM* is the best constrained quantity for the majority of 
objects; we refer to the above papers for a more detailed dis- 
cussion 13 . The derived uncertainties are smaller or compa- 
rable to the differences in best-fit values obtained by va rying 
the assumed SFH, metallicity, or extinction law (see § 15. 4> . 
The main goal of this work is to characterize the ensemble 
properties of DRGs (from the median and the global distribu- 
tions of best-fit values), and accounting for the uncertainties 
of individual sources does not alter our general conclusions. 
We will therefore not discuss them further, except where most 
relevant. 

5.3. Modeling Results for the Standard Model Sets 
Both the CSF and T3<joMyr models provide acceptable fits to 
the SEDs of our DRG samples 14 . The best-fitting models for 
selected sources are shown in Figure E| Tableland Fi gureiyj 
summarize the derived properties. The results are similar for 
the HDF-S and MS 1054-03 field samples. With the CSF 
model, the median values for HDF-S (MS 1054-03 field) are 
1.7 (2.0) Gyr, (Ay) med = 2.7 (2.4) mag, (M*) med = 
8.1 (16) x 10 10 M , (M+fLv^med = 1-2 (2.3) MqL^q, and 
instantaneous (SFR) med = 120 (170) M Q yr _1 . With the 

13 For the combined DRG sample, the median of the uncertainties given 
by the 68% confidence intervals for the CSF <T3ooMyr) standard model are 

^% (+|8%) for the age frf , +26% (+28%) for M ^ +14% ( +2l%) for M*/L y ,*, and 

^j>* (50% to a factor of 2) for the instantaneous SFR. For the HDF-S LBG 
sample also modeled in this work ($ 16.11 . the uncertainties are C29*) 
for the age t a , + 1™ C' 3 *) for M*, +22% (±5%) for M„/L v ^ and « 3 « 
for the instantaneous SFR. For DRGs and LBGs, the uncertainties on Ay are 
< 0.3 mag with both the CSF and T3ooMyr models. 

14 The median of the chi-squared values normalized per degree of freedom, 
Xn. for the best-fit CSF models are 1.6 for the HDF-S sample and 2.1 for the 
MS 1054—03 field sample. For the T3ooMyr models, the values are 0.6 and 
1.6, respectively. Xn' s higher than the expected value of unity for good fits 
could indicate that the models do not provide a realistic representation of the 
objects, the photometric uncertainties underestimate the real measurement 
errors, the distributions of uncertainties of the modeled properties are not 
Gaussian, or, most likely, that a combination of the factors is at play. 



T3ooMyi model, the median age and extinction are 1.0 Gyr 
and 0.9 mag for both fields; the median for the other prop- 
erties in HDF-S (MS 1054-03 field) are 7.5 (10) x 10 10 M Q , 
1.2 (1.8) M Q Lv| Q , and 7 (23) M Q yr _1 . The median ages of 
the bulk of stars (f s f)sFR and the median time-averaged SFRs 
M*/f s f for the two SFHs and two fields differ less than the 
median f s f and instantaneous SFRs, with values between 0.75 
and 1.0 Gyr, and 90 and 140 Moyr" 1 . 

The dispersion among individual DRGs is large for all prop- 
erties. This stems partly from the diversity in SEDs noted in 
§ |4] For instance, with a given SFH, quite different com- 
binations of age and extinction are derived for the objects 
that are bluest in the rest-frame UV compared to the red- 
dest ones. Part of the (small) differences between the HDF- 
S and MS 1054-03 samples likely reflects the brighter K s - 
band magnitude limit and larger proportion of the brightest 
sources, the redder J s -K s colours, and the lower redshifts for 
the MS 1054-03 objects (see Table [l). The different sets of 
optical photometric bands may play a role too by providing 
somewhat different observational constraints. 

For the ensemble of DRGs, either SFH thus implies old 
ages, and high stellar masses and M/L ratios. High SFRs and 
large amounts of extinction are derived for the CSF model 
while the T3ooMyi model leads to lower SFRs and extinction 
(we c omment on the validity of the adopted ICalzetti et al. I 
2000 extinction la w in v iew of the large extinction derived 
for the DRGs in § 15.4.21 . Simila r results were obtained in 
our initial study of HDF-S ( Franx et al~l 120031) and for the 
spectroscopic sample (van Dokkun i et al. 1 12004 > . In particu- 
lar, for the two spectroscopically-confirmed non-AGN DRGs 
in the MS 1054 -03 field, the ages are supported by esti- 
mates from the Ha equivalent widths and the stellar masses 
are consistent with the dynamical masses from the Ha or 
[O III] A4959,5007 A line widths. As we show in § EU the 
stellar masses are the most robust properties against variations 
in SFH, metallicity, and extinction law (neglecting changes in 
the IMF). With our definition, the M*/Ly* ratios also are be- 
cause the luminosities involved are extrinsic. On the other 
hand, the instantaneous SFRs are very sensitive to the model 
assumptions, especially to the SFH (again ignoring the IMF). 

To test the sensitivity to the photometric redshift uncertain- 
ties, we modeled the DRG samples as above, adopting in turn 
the lowest and highest redshift allowed by the 68% confidence 
interval around z p h for each object. Although the effects on 
the median values of the best-fit properties can be appreciable 
(by up to 0.7 Gyr for the age, 0.9 mag for Ay, a factor of 2 
for M+, and a factor of 5 for SFR, depending on the field and 
SFH), our conclusions are qualitatively unaffected. We also 
modeled the samples of DRGs including the objects outside 
of the z = 2-3.5 range. The impact on the median properties 
is significantly smaller (< 0.5 Gyr, < 0.3 mag, < 15%, and 
< 30% for the ages, Ay's, M*'s, and SFRs, respectively). 

The SED fits of Figure [4] illustrate well the difficulty in 
discriminating between the two SFHs (in particular, e.g., for 
HDFS-496 and MS - 1383). Based on the formal x 2 val- 
ues, the fits are statistically indistinguishable for all but three 
galaxies, for which the CSF model could in principle be ruled 
out at the 95% confidence level. For two of these, HDFS-767 
and MS - 1038, the higher flux in the H band relative to the 
K s band strongly influences this result. Although their H-K s 
colour is fairly blue, they are not the most extreme of our 
sample. For the other, MS - 1319, the optical data are bet- 
ter reproduced with the T3ooMyr model. For this spectroscop- 
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ically studied DRG ( Ivan Dokkum et al. 12004) . correcting for 
line emission would affect mostly the observed Z?-band flux 
(downwards by 10% for Lya) and would improve the agree- 
ment with the CSF model. The distinction between the mod- 
els is not very robust and the \ 2 computation relies on uncer- 
tainties that do not account for those of z p h or possible emis- 
sion line contributions. In view of this, we chose not to dis- 
criminate between SFH models for individual sources in our 
analysis. 

5.4. Variations in Model Parameters 

We explored the effects of variations in three of the model 
assumptions: the SFH, the extinction law, and the metallic - 
ity. Figure^] summarizes the results for the ages, extinction, 
stellar masses, and SFRs derived for different sets of input 
parameters. While the effects may differ among individual 
objects, we are here interested in the ensemble properties of 
DRGs and so the discussion refers to the median values. We 
find that the SFH has the largest impact on the derived quanti- 
ties. The age, extinction, and instantaneous SFR are the most 
affected ones whereas the stellar mass and M/L ratio are the 
least sensitive ones and change only by factors of m 2-3. 

5.4. 1 . Star Formation History 

We considered exponentially declining SFRs with a range 
of e-folding timescales r between 10 Myr and 1 Gyr. They are 
compared with the CSF and also an SSP model, equivalent to 
r = oo and r = 0, respectively. With our broad-band data, it 
is not possible to discriminate among any of the SFHs for the 
large majority of our DRGs. The overall trend with decreasing 
timescale is of younger ages by up to an order of magnitude. 
The extinction values become lower but do not vary mono- 
tonically with r, exhibiting a minimum around 300 Myr. The 
maximum difference amounts to 1.8 mag, or a factor of « 5 
in attenuation at V in the rest-frame. It is intriguing that the 
lowest extinction values are not derived for the models with 
shortest timescales, as might be expected. This may be due 
to worse template mismatch at r < 10 8 yr. The SFRs drop by 
several orders of magnitude, to some extent because of lower 
extinction but mainly as a consequence of the functional form 
of the assumed SFR itself combined with the ages, which are 
such that R(t s f)/Ro = exp(-f sf /r) < 1. 

The stellar masses tend to decrease with smaller r but re- 
main within a factor of 2 of the masses for CSF. We thus 
reach similar concl usions about th e stability of deriv ed M f 
for DRGs as Pa povich et ail d2001l) and lShaplev et al. I d200lfe 
in their modeling of LBGs also based on broad-band optical- 
to-NIR photometry. As discussed by these authors, this is 
in part due to the NIR data probing the rest-frame optical 
regime, which is less sensitive to M/L variations due to age 
or extinction than the rest-frame UV. In addition, the age and 
dust degeneracies tend to cancel out in terms of the effect 
on the ratio of stellar mass to observed (attenuated) light. 
The same behaviour is seen in the case of our DRGs (see 
Ivan Dokkum et ain i2004. for a discussion). While age and 
extinction can vary strongly with the SFH, the combinations 
of best-fit values result in much tighter constraints on and 
M*/L v .*. 

The largest extinction is derived for the longest timescales. 
The need for high extinction here may be artificially enhanced 
by the choice of SFH along with the constraint that the best- 
fit ages do not exceed that of the universe at each object's 
redshift. In the limit of CSF, the population of young mas- 
sive stars rapidly levels off on the timescale of their main- 



sequence lifetime (~ 10 6 - 10 7 yr) while the population of 
lower-mass longer-lived stars grows progressively. To match 
the typically red SEDs of DRGs, the UV light from the lumi- 
nous massive stars continuously formed has to be suppressed 
with sufficient extinction. The maximum age constraint lim- 
its the relative contribution of evolved versus young stars to 
the integrated light, further enhancing the need for extinction. 
Over the range of SFHs explored here, the strong coupling 
between age and r also reflects the need for the build-up of a 
large enough population of evolved stars. 

Models with very short timescales < 10 8 yr where the ex- 
tinction is allowed to vary lead to ages < 300 Myr and fairly 
high AyRi 1.5 mag. Arguably, such young dusty burst mod- 
els may not be appropriate for the ensemble of DRGs (though 
they may be for some individual objects) in view of the evi- 
dence for a prominent spectral break in their rest-frame con- 
tinuum SEDs (see §|4j. Choices of r ~ 10 s yr lead to ages 
consistent with the appearance of a strong Balmer/4000A 
break in both cases of very quiescent or still actively star- 
forming galaxies. Although models with t < 10 8 yr cannot 
be statistically ruled out in our fits, we regard them as least 
plausible for DRGs. 

5.4.2. Extinction Law and Metallicity 

To test the effects of possible variations in the extinction 
law and in metallicity, we fitted two alternative suites of mod- 
els to the SEDs. In one variant, we modifie d the p rescrip- 
tion for the extinction law to that of Allen ( 1976) for the 
Milky Way (MW), keeping the solar metallicity. In the other, 
we used the set of BC03 models for metallicity Z = 0.2 Z Q 
and adopted for consisten cy the extinction law for the Small 
Mage llanic Cloud (SMC; iPrevot et al~lll984t iBouchet etai~| 
1985). The main differences between the ICalzetti et al. I 
(2000) and MW extinction laws lie in the ratio of total-to- 
selective absorption Ry = Ay /E(B -V) (4.05 versus 3.1, re- 
spectively) and in the ICalzetti et al. I law lacking the 2175 A 
bump characteristic of MW dust mixtures. Otherwise, their 
wavelength dependence are fairly similar. The SMC law with 
Rv = 2.72 also lacks the 2175 A bump. In addition, it rises 
more steeply with decreasing wavelengths i n the near-UV 
than the other two laws; in other words, the ICalzetti et al. I 
and MW laws are much "greyer" at near-UV wavelengths. 

For the models with Z = Z Q and the MW law, the ages for 
a given SFH remain the same or become mostly younger, by 
up to a factor of 3. The best-fit Ay's also remain the same or 
decrease, by at most 0.6 mag. The effects are more important 
for the Z = 0.2 Z Q and SMC law models, with generally older 
ages by up to a factor of 5 and Ay's lower by up to 2 mag. 
In both model sets, the instantaneous SFRs are strongly af- 
fected towards shorter timescales t, especially for the low- 
metallicity case. The stellar masses are systematically lower 
for both variants, by a factor of « 2 on average. 

The near-systematic increase in median age of the DRGs 
for the lower metallicity models can be attributed to the com- 
bined effects from the extinction law and stellar emission. 
Si nce the SMC law is ste eper in the rest-frame UV than both 
the Calzetti et al.l J2000) and MW laws, lower Ay values will 
reproduce the observed SED slopes in this regime but also 
imply bluer rest-frame UV-to-optical colours such that older 
ages are required to match the colours of DRGs. In addi- 
tion, at Z < Z Q and for a fixed age, the stellar evolutionary 
tracks and model atmospheres predict higher effective tem- 
peratures, bluer continua, and shallower metallic and molec- 
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ular absorption features, furth er driving up the b est-fit model 
ages. As emphasized by Shaplev et al J J2001I) . it is impor- 
tant here to keep in mind that theoretical stellar tracks and 
stellar atmosphere models for non-solar metallicities are still 
not fully tested against e mpirical data (see the discussion by 
Bruzual & Chariot 2003). The variations in derived parame- 
ters we describe above provide nevertheless indications of the 
magnitude of metallicity effects. 

We further note that the high obscuration inferred for the 
models with star formation timescales r > 10 9 yr raises 
concerns abou t the applicability to DRGs of the adopted 
iCalzetti et al. I (120001) extinction la w and foreground dust 
screen geometry. The ICalzetti et al. I law provides a good rep- 
resentation for local UV-bright starbursts with Ay up to about 
2.5 mag and where the foreground dust approximation is ob- 
servationally supported. For the CSF model, for instance, the 
Ay values obtained for half of the DRGs indicate extinction 
levels at the limit of , or beyond the range of validity. This 
may suggest that the Calzetti et al. law is inappropriate, and 
that the dust and sources geometry is more complex. 

6. COMPARISON OF DRGS WITH LBGS AND OTHER 2 < Z < 3.5 
GALAXIES 

The picture that emerges from our SED analysis and mod- 
eling is that DRGs constitute an old and massive galaxy pop- 
ulation at high-redshift. Obscuration by interstellar dust plays 
an important role in their observed properties. Depending on 
their actual SFH, they may also still be very actively forming 
stars. In this section, we interpret our results in a broader con- 
text by comparing our DRGs with LBGs. We also consider 
the ensemble of all /Q-band selected galaxies in both FIRES 
fields at 2 < z < 3.5. We modeled these samples in the exact 
same way as the DRGs. 

6.1. How Different are DRGs and LBGs ? 

An immediate implication of the model results presented in 
§|5]is that DRGs appear to be older, more obscured, and with 
higher stellar masses than LBGs. Our sample of optically- 
selected LBGs at Kl ot < 22.5 mag and 2 < z < 3.5 in HDF- 
S allows a consistent comparison, based on the same data 
and with identical model ingredients and assumptions. With 
our standard CSF model, we derive median values for the 
LBGs of (f s f} me d = 500 Myr, (A v ) med = 0.6 mag, (M*) med = 
1.4 x 10 10 M Q , {M,/L v .+) med = 0.4 M Q L V ; Q , and (SFR) med = 
39 M Q yr _1 . With the T3ooMyr model, we find lower median 
age of 360 Myr and SFR of 21 M© yr" 1 , but the same median 
extinction, stellar mass, and M/L ratio. We cannot statistically 
discriminate between the fits for different SFHs for any of our 
LBGs. 

Compared to the HDF-S LBGs and with CSF, our DRGs 
are typically 3-4 times older, 4-5 times more obscured, have 
roughly an order of magnitude higher stellar mass, and form 
stars at 3 - 5 times higher instantaneous rates (based on the 
median properties). Among the various declining SFHs we 
considered, the T3ooMyr model leads to the lowest extinction 
for the DRGs (Figure [Toj. One might expect that this SFH 
reduces most the inferred differences between the two pop- 
ulations. The DRGs, however, remain distinct from LBGs 
in their median properties except for the SFR, being 3 times 
older, 1.5 times more extincted, 5-8 times more massive, and 
0.3- 1.3 times as actively star-forming. As for the DRGs, the 
LBGs show an appreciable dispersion in best-fit properties 
with some overlap with DRGs. However, the differences in 
the median values and combinations thereof (especially age, 



Ay, and M*) indicate that the ensemble properties of the two 
populations are different for either SFH. 

The distinction between the DRGs and LBGs is best illus- 
trated in Figure which compares the distribution of the 
two populations in the M^/Ly^-Ay plane (where the median 
uncertainties from the 68% confidence intervals are < 25% in 
M*/Ly,* and < 0.3 mag in Ay for the DRG and LBG sam- 
ples, and our two standard models). In this parameter space, 
the overlap is very small for the CSF model. For the T3ooMyr 
model, about two-thirds of the DRGs shift to a lower ex- 
tinction interval Ay < 1.5 mag, in common with the LBGs, 
while the rest uniformly populates the higher Ay range ex- 
plored. The ensemble of DRGs and LBGs remain however 
clearly separated in the ranges of M+/Ly^ ratios covered, 
with very little overlap. This is driven by the stellar masses 
since the DRG and LBG samples are similarly luminous in 
the rest-frame V band, with nearly identical median values of 
5 x 10 10 L v ,o (within 20%). 

The weakness of our HDF-S LBG sample is that z sp 's are 
available for seven out of 33 galaxies. However, we arrive at 
similar conclusions if we consider the model results for the 
spectro sco pically confirmed LB G samples of Pa povich et alJ 
( 2001 ) and Shap lev et al. H2001I) . also based on optical to NIR 
SEDs. The typical (median or geometric mean) ages inferred 
are ~ 10 s yr, with extinction Ay < 1 mag andM* ~ 10 10 M Q 
(for Z = Z , a Salpeter IMF down to 0. 1 M©, and our adopted 
co smology). Our r esul ts are within a f actor of a few of those 
of Pa povich et alJ and Shaple vet al. L similar to the differ- 
ence between the latter two. It is however difficult to com- 
pare in detail the results among all three studies. The data 
sets are different in terms of depth, quality, and photomet- 
ric bands involved. The sample of Papovich et al. (2001) 
was selected in the HDF-N from the space-based HST data 
while the ground-based sample of Sha plev et al. I ll2001l) was 
select ed over a much wider area and is brighter. In their mod- 
eling,|Papovich et al. (2001) assumed exponen tially declining 
SFHs with the timescale as a free parameter. Shaple vet al. I 
(2001) did not attempt to fit the SFH and discuss mainly 
the CSF scenario. In addition, the models of Papovic h et alJ 
and Shaple vet al. I relied on earlier versions of the Bruzual & 
Chariot code and partly different model ingredients. In view 
of this, it is probably fair to say that the properties derived for 
all three LBG samples are in broad agreement. The conclu- 
sions about the differences between DRGs and LBGs remain 
qualitatively the same. 

6.2. Derived Properties versus Observed NIR Properties 

We investigated relationships between the modeled quanti- 
ties and the observed NIR properties, the key feature of our 
DRG selection. For this purpose, we complemented the DRG 
samples with the HDF-S LBGs and all other /^ s -band selected 
galaxies at 2 < z < 3.5 in both the HDF-S and MS 1054-03 
fields down to the same limits of 22.5 and 21.7 mag, re- 
spectively. 

Figure El plots, for our standard T3ooMyr model, the de- 
rived stellar masses and M/L ratios as a function of observed 
J s -K s colour and X" s -band magnitude. The results with the 
CSF are qualitatively similar. There is a clear correlation be- 
tween the derived and the observed J s -K s colour, which is 
even tighter for the M+/Ly^ ratio. The redder the objects, the 
higher the stellar masses and mass-to-light ratios. The DRGs 
lie at the massive end and extend smoothly the relationships 
for the 2 < z < 3.5 objects with bluer J s -K s colours, includ- 
ing the LBGs. Because of the different limiting magnitudes 
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between the HDF-S and MS 1054-03 fields, it is difficult to 
draw conclusions on the variations of and M+fLv.* ver- 
sus Kl° l . The data suggest that our ,ff s -band selection does not 
miss a significant fraction of massive galaxies at 2 < z < 3.5 
that would be faint because of extinction effects. 

Figure^]plots the variations of the derived extinctions and 
instantaneous star formation rates per unit stellar mass with 
observed J s -K s colours. Results for both the standard CSF 
and T3Q0Myr models are shown. There are general trends of in- 
creasing Ay and of decreasing SFR/M* with redder J s -K s 
colour. The variations in derived ages (not shown in Fig- 
ure are essentially the reverse of those in the SFR/M* 
ratio, which provides a measure of l/f s f (with a direct pro- 
portionality for the CSF model, and a non-linear monotonic 
relation for the T3ooMyr model). These plots further show that 
if DRGs might be exceptional in their high absolute SFRs, 
they are not at all in terms of their SFR per unit mass. 

6.3. Past History and Future Fate of DRGs 

We here speculate about the evolution of the 2 < z < 3.5 
DRGs in the FIRES fields, based on the results for the CSF 
and T3ooMyr models. Figure plots the relative variations 
with redshift of the average SFR of the DRG samples and, for 
comparison, of the HDF-S LBG sample at 2 < z < 3.5. We 
computed the curves from the model SFR R(t) used as input 
in the evolutionary synthesis. For each galaxy, we scaled the 
model R(t) curves to match the derived instantaneous SFR at 
the epoch of observation. For the purpose of this discussion, 
we in fact assumed that all galaxies lie at z = 2.5, close to the 
mean and median redshift of both our DRG and LBG sam- 
ples (see §0. We then normalized the resulting average SFR 
curves to unity at z = 2.5. 

For the CSF model, the average instantaneous SFR of the 
DRG samples has not varied strongly since z = 6 up to z = 2, 
increasing by a factor of ~ 2 with time. On the other hand, the 
average instantaneous SFR of the HDF-S LBG sample has in- 
creased by more than an order of magnitude over the same 
redshift interval. For the T3ooMyr model, the variations are 
much noisier but indicate, broadly speaking, roughly com- 
parable activity levels in the past out to z « 4, where the 
mean instantaneous SFR of LBGs drops abruptly. The rel- 
ative variations in average SFR backwards in time are domi- 
nated by the age distribution among each population, and the 
steeper drop for the LBGs reflects their younger ages com- 
pared to the DRGs. This extrapolation of the star formation 
histories suggests that 2 < z < 3.5 DRGs have started form- 
ing stars at very early epochs, and earlier than 2 < z < 3.5 
LBGs have. A possible scenario for DRGs is that they 
are the evolved descendants of systems that experienced a 
"Lyman-break phase" at z > 4 (see also iFranx et al. 1 120031: 
Ivan Dokkum et akl 120031 120041) . This evolutionary picture 
may be consistent with the increased redde ning in rest-fr ame 
UV- to-optical colours inferred by Pacovichet al. (2004J) for 
Lyman-break selected galaxies between z ~ 4 ("Z?-dropouts") 
and z ~ 3 ("{/-dropouts"). We emphasize that the above dis- 
cussion applies to the stellar populations that dominate the 
rest-frame UV-optical emission of the galaxies, as traced by 
the observed SEDs. 

As discussed by van Dokkum TTaH (|2004), the stellar 
masses of DRGs are already very high atz ~ 2-3 and compa- 
rable to those of early-type galaxies locally. This may indicate 
that DRGs have already assembled most of their stellar mass 
at these redshifts and evolve thereafter mostly passively. Us- 
ing our CSF model results, we estimated the stellar masses 



predicted at z = if star formation is maintained at the same 
rate as inferred at the redshift of observation. The median 
present-day stellar masses would be 2 x 10 12 M Q , which cor- 
respond to the most massive galaxies in local galaxy clusters. 
While some of the DRGs might evolve into such systems, it 
seems rather unlikely we have detected such large numbers 
of their progenitors in our two small fields. More plausibly, 
most of the DRGs probably become predominantly quiescent 
systems at z < 2, and may possibly grow further in mass dur- 
ing brief merger events. Assuming that DRGs stop forming 
stars around their redshift of observation and evolve passively 
thereafter, the CSF models predict colours of I-K s > 4 and 
R— K s > 5 at 1 < z < 2, even if all the dust is removed at the 
instant when star formation is quenched (such red colours are 
also predicted using the T3ooMyr models instead). This may 
suggest an evolutionary connection between DRGs at z > 2 
andEROs atz < 1.5. 

In absolute terms, the instantaneous SFRs of ~ 100 M Q yr" 1 
derived from the SED modeling assuming CSF lie at the high 
end inferred for the local galaxy population, in the regime 
of luminous and ultra-luminous infrared galaxies. Such a 
high SFR rate seems supported by the possible detection of 
one of the non-AGN spectroscopically-confirmed DRG in the 
MS 1054-03 field at sub-millimet er wavelength s from deep 
SCUBA observations ( Ivan Dokkum et al. 1120041 K. Kraiberg 
Knudsen et al. 2004, in preparation). The CSF hypothesis is 
unverifiable with the data at hand, and may be a valid repre- 
sentation for a few objects only. Conversely, while the T3ooMyi 
model implies substantially lower instantaneous SFR at the 
redshift of observation, the high stellar masses lead to ex- 
tremely high SFRs at earlier epochs (> 1000 M Q yr -1 ). Obvi- 
ously, the SFH of DRGs may well be more complex than our 
simplistic assumptions. The less model-dependent estimate 
fromM^/fsf of ~ 100 M yr _1 seems to suggest elevated lev- 
els of star formation activity averaged over the lifetime of the 
stellar populations traced by our SEDs. 

6.4. Relative Importance of DRGs and LBGs at 2 < z < 3.5 

Using the results presented in this paper, we can crudely 
assess the relative importance of the DRG and LBG popula- 
tions at high redshift. In the following, we consider our sam- 
ples of DRGs, LBGs, and Zf s -band selected sources lying at 
2 < z < 3.5 and to the K s total magnitude limits of 22.5 for 
HDF-S and 21.7 for the MS 1054-03 field. In HDF-S, the 
LBGs outnumber the DRGs by a factor of three. They con- 
tribute about 2.5 times more than DRGs to the total observed 
/iT-band light and extrinsic rest-frame V-band luminosity Ly,* 
from all sources, with fractions of sa 60% and « 25%, respec- 
tively. The relative contributions are reversed for the stellar 
mass, where DRGs account for w 60% of the total compared 
to w 35% for LBGs (using the individual M + 's derived from 
either the CSF or T3ooMyr models, and assuming the same 
model parameters for all populations). In the MS 1054-03 
field and to the corresponding brighter magnitude limit, the 
DRGs produce about 35% of the total observed /f s -band light 
and Lv,±, and make up m 65% of the integrated M*. 

This provides a first-order indication that DRGs may repre- 
sent a significant constituent of the z ~ 2 - 3 universe in terms 
of stellar mass. Admittedly, our estimates are rough. A de- 
tailed analysis of the contribution of DRGs to the stellar mass 
budget in both FIRES fields wi ll be presented by G. Rudnick 
et al. (in preparation; see also Rudnick et al. 2003 for first 
estimates in HDF-S). 
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7. SUMMARY AND CONCLUSIONS 

We have presented the analysis and theoretical modeling of 
the A = 0.3-2.2 ptm SEDs of 34 DRGs at 2 < z < 3.5. The 
galaxies were selected by their J s -K s > 2.3 colours in the 
HDF-S and MS 1054-03 fi elds observed as part of the FIRES 
survey jFranxet al.ll20 00). at K s < 22.5 and < 21.7 mag, re- 
spectively. The addition of the MS 1054-03 field increases 
the area surveyed by a factor of five to an « 0.7 mag brighter 
limit and triples the original sample from HDF-S. This al- 
lows for the first time a robust assessment of the ensemble 
properties of DRGs. We find large numbers of J s -K s red ob- 
jects in both fields, with surface densities of ~ 1 arcmin -2 at 
K s < 21 mag. Since no J s -K s > 2.3 source is found in the 
HDF-N down to K s = 21 mag, cosmic variance is substantial. 
Observations of wider areas and multiple lines of sight will 
be required for a more reliable determination of the surface 
density of the J s -K s red population. 

Using new diagnostic diagrams involving the 1^4 -J s , J s — 
H, and H - K s colours, we have shown that the red NIR 
colours of the DRGs cannot be attributed to extinction effects 
alone and require the presence of a prominent Balmer/4000 A 
break characteristic of evolved stars for a large fraction of 
them. In the rest-frame optical, the SEDs of DRGs fall 
within the envelope of normal nearby galaxies and, in the 
rest-frame UV, indicate a wide range in star formation activity 
and/or dust obscuration. The rest-frame UV-to-optical SEDs 
of DRGs are in stark contrast with those of LBGs, which are 
much bluer and more uniform. These differences can be eas- 
ily understood in terms of selection effects. The LBG crite- 
ria preferentially select galaxies that are actively star-forming 
and modestly obscured. Our J s -K s > 2.3 criterion primarily 
targets evolved systems, ranging from those with quiescent or 
no star formation and little extinction to those that are highly 
obscured and vigorously forming stars. 

We have applied evolutionary synthesis models to the SEDs 
of our DRGs to constrain quantitatively the age of the stellar 
populations, the extinction, the stellar masses, and the star 
formation rates. We have also investigated the effects of vari- 
ations in the assumed star formation history, metallicity, and 
extinction law on the model results. The derived properties 
are similar for the HDF-S and MS 1054-03 field samples. 
For constant star formation, a Salpeter IMF betwe en 0.1 and 
100 M Q , solar metallicity, and the lCalzetti et al. I J20QOI) ex- 
tinction law, we obtain for each field a median age ~ 2.0 Gyr, 
extinction Ay ~ 2.5 mag, stellar mass ~ 10 11 M©, stellar 
mass-to-light ratio M+/Lv.+ ~ 1 M Q Ly 1 Q , and instantaneous 
star formation rate SFR ~ 100 Moyr" 1 . The stellar masses 
and M/L ratios are the most robust properties against varia- 
tions in the model parameters (neglecting those of the IMF) 
while the instantaneous SFRs are strongly model-dependent. 
Models with exponentially declining star formation rates and 
e-folding timescales in the range r = 10 Myr- 1 Gyr generally 
lead to lower ages, extinction values, and instantaneous SFRs 
but similar stellar masses and M/L ratios within a factor of 
two. The smallest median Ay's are derived with r = 300 Myr. 
For this T3ooMyr model and the assumptions above for the 
IMF, metallicity, and extinction law, the median age and Ay 
are 1 Gyr and 0.9 mag for both fields, and the median SFRs 
^lOMoyr" 1 . 

An important implication of our work is that DRGs appear 
to be dominated by significantly older, more obscured, and 
more massive stellar populations than LBGs at similar red- 
shifts and similar rest-frame V-band luminosities. Using our 



sample of optically-selected LBGs in HDF-S, we have ex- 
plicitely shown that the distinction holds for constant star for- 
mation and T3ooMyr models, the two star formation histories 
that bracket the range of median Ay's derived for our DRGs 
(neglecting variations in model parameters between the two 
classes). For both DRGs and LBGs, there is a large spread in 
derived properties among individual objects and some of the 
DRGs have ages, extinction, and stellar masses characteristic 
of LBGs, and vice-versa. However, the overlap is small, espe- 
cially for the stellar masses and M/L ratio s. C omparison with 
the LB G studies of Papo vich et alJ J2001I) and Shaple vet al. I 
( 2001 ) is less straightforward because of differences in the op- 
tical/NIR data sets, in the details of the sample selection, and 
in the modeling but, broadly speaking, it supports the same 
conclusions. More generally, we find that the derived age, ex- 
tinction, and stellar mass all increase with redder J s -K s colour 
among the full sample of 2 < z < 3.5 K s -band-selected galax- 
ies in the FIRES fields, down to the same magnitude limits 
as the DRG samples. The relationships are particularly well- 
defined for the stellar mass and M/L ratio versus J s -K s colour. 

If DRGs represent the most evolved and most massive 
galaxies at z ~ 2.5, it is tempting to speculate that they might 
be the precursors of the most massive present-day galaxies. 
Such a final fate has been proposed for LBGs at z ~ 3 in 
the context of hier archical models of galaxy formation (e.g. 
Baugh et al. 1998). For DRGs, with stellar masses already 
similar to those of local early-type galaxies, one may rea- 
son in a simplistic way that there is little room for substan- 
tial growth through subsequent merger events. Extrapolating 
the star formation histories backward in time for both DRGs 
and LBGs at 2 < z < 3.5, we infer that the DRG population 
has been forming stars at a roughly constant rate and possibly 
back to z ~ 5 - 6, whereas the LBGs started to form later at 
z ~ 4. This is consistent with a scenario in which the DRGs 
are the evolved descendants of systems undergoing a "Lyman- 
break phase" at z > 4. The DRGs may evolve into the ERO 
population at 1 < z < 2. Another intriguing possibility is a 
connection with the high-redshift population detected at sub- 
millimeter wavelengths. One of the DRGs in the MS 1054-03 
field likely is the cou nterpart of a sub-m illimeter source de- 
tected with SCUBA llvan Dokkum et al. II200H K. Kraiberg 
Knudsen et al. 2004, in preparation). In addition, although de- 
tailed information is currently available for a few objects only, 
recent work indicates that some of the bright sub-millimeter 
galaxies at z = 2 - 3 host evolved stellar populations, show 
evidence of substantial chemical enrich ment, and have dy- 
namical/stellar masses ~ 10 11 M^ (e.g. iGenzel et al. Il2003t 
iNeri et al. 1200 31 iTecza et al.l200l . 

The results presented here constitute a major improvement 
over our previous studies of J s -K s selected objects at z > 2. 
Significant uncertainties remain, however, and further inves- 
tigations are clearly needed. Our analysis still relies on a lim- 
ited sample selected from two small fields and the model pa- 
rameters are largely unconstrained. Deep photometric surveys 
of wider areas and disjoint fields as well as extensive follow- 
up spectroscopy will be crucial to assess (and average over) 
field-to-field variations, secure the redshifts, derive robust 
number densities, and constrain more accurately the physical 
properties (e.g., extinction, metallicity, dynamical masses). 
The current dearth of constraints on the nature of DRGs fur- 
ther limits our ability to interpret their observed properties. In 
that respect, it will be very important to establish firmly the 
frequency of AGN among DRGs and their contribution to the 
measured optical/NIR SEDs, in particular from spectroscopy. 
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It will also be important to obtain independent constraints on 
the current and past star formation rates, with observations at 
infrared and submillimeter wavelengths (e.g., with SCUBA 
and the Spitzer Space Telescope) complementing rest-frame 
optical line diagnostics. Ultimately, a better understanding 
of DRGs and their relationship to other galaxy populations at 
high and low redshift may provide important constraints on 
the assembly history of massive galaxies. 
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FIG. 1. — Redshift distributions of the J s —K s > 2.3 selected galaxies in the FIRES fields. The sample includes objects at A^ ot < 22.5 mag from HDF-S 
and K'" 1 < 21.7 mag from the MS 1054 — 03 field, (a) Comparison between the photometric and spectroscopic redshifts for the objects with z sp determinations 
(available in MS 1054-03 only). The error bars represent the Icr uncertainties on z p h- The diagonal line indicates direct proportionality between z f h an d Zsp. 
Object MS — 140 is also show n he re (open circle) although it is not formally part of the sample; it satisfies all the selection criteria except the minimum exposure 
time in the NIR bands (see § 13.21 . (b) Redshift distributions of the HDF-S (hatched histogram) and MS 1054 — 03 (grey-filled histogram, excluding MS - 140) 
samples. When available, the spectroscopic redshift was used instead of the photometric redshift. 
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FIG. 2. — Cumulative A^-band raw number counts of / s — K s > 2.3 mag objects in the FIRES fields. The counts for HDF-S (triangles) and the MS1054 — 03 
field (circles) include J s —K s selected objects down to the respective 5<r total limiting magnitudes for point sources of A"J ot = 23.8 and 23.1 mag. The effective 
selection areas are 4.48 arcmin 2 for HDF-S and 23.86 arcmin 2 for the MS 1054 — 03 field (see § 13.21 . The counts for the MS1054 — 03 field are not corrected for 
the lensing effects of the cluster, which are small (on average ?s 20% in both flux and area magnification) and to first order cancel out. Filled symbols indicate 
counts for objects with redshift z, > 2 while open symbols show counts with no redshift restriction (the z sp was used whenever available). The error bars indicate 
the Poisson uncertainties. The vertical lines show the magnitude cutoffs of the samples adopted for analysis ensuring a S/N > 15 on the S" s -band colour flux, at 
Xf' = 22.5 for HDF-S (dotted line) and 21.7 mag for MS 1054-03 (dashed line). 
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FIG. 3. — Cumulative raw number counts as a function of Js—Ka colour cutoff for Xs-band selected sources in the FIRES fields. The counts are given as number 
of objects per unit area N with colour redder than the J s —K s of the abscissa (the counts for the MS 1054-03 field are uncorrected for the area magnification 
by the cluster, a small effect of a 20% on average). The counts for HDF-S (triangles) and the MS 1054 -03 field (circles) include objects down to A^" 1 = 22.5 
and 21.7 mag, and over the effective selection areas of 4.48 and 23.86 arcmin 2 , respectively (see § 13.21 . Filled symbols indicate counts for objects with redshift 
2 < z < 3.5 while open symbols show counts with no redshift restriction (the z sp was used whenever available). The error bars indicate the Poisson uncertainties. 
The vertical dotted line shows the J s —K s colour cutoff at 2.3 mag. Most galaxies with J s —K s > 2.3 colours have redshifts at z > 2. 
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FIG. 4. — Broad-band SEDs of selected / s — ^ s > 2.3 objects at 2 < z < 3.5 in the FIRES fields. The panels to the left show SEDs of HDF-S objects and to the 
right, SEDs of objects in the MS 1054—03 field. The galaxies are sorted by increasing /gi4 -K s colour. All SEDs are corrected for Galactic extinction. The ve rtica l 
error bars and the upper limits represent the lcr uncertainties on the flux measurements or the minimum error of 0.05 mag assumed for the modeling (see t? !5 . 21 . 
The h orizontal error bars indicate the bandpass widths at half the maximum transmission. The best model fits for the two "standard" star formation histories (see 
S I5.1I are overplotted: a constant star formation rate (CSF; black lines) and an exponentially declining star formation rate with r = 300 Myr (T3ooMyr; 8 re y lines). 
The best-fitting age and visual extinction Ay are given for both models in each panel. The models were computed for a Salpeter IMF between 0.1 and 100 Mq, 
solar metallicity, and the Calzetti et al. 1 2000) extinction law. 
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FIG. 5.— Rest-frame broad-band SEDs of 2 < z < 3.5 galaxies in the FIRES fields, (a) J s -K s > 2.3 selected objects at Kl°' < 22.5 mag in HDF-S 
and Kl ot < 21.7 mag in the MS 1054 — 03 field. The objects with available spectroscopic redshift are plotted with large white-filled dots (including MS — 140, 
see § 13.21 . (b) Optically-selected LBGs at K ^ ot < 22.5 mag in HDF-S. All SEDs are corrected for Galactic extinction, and for the average absorption due to 
intergalactic Lyman line blanketing following Madau 1 1995). The error bars represent the lcr uncertainties on the flux measurements while upper limits indicate 
the 2cr confidence levels. The solid lines show the template spectra for galaxy types E, Sbc, Scd, and Im from Coleman et al. 1 1980), extended beyond the 
original wavelength range using evolutionary synthesis models. The shaded area outlines the range in SED shapes covered by these templates within the original 
wavelength limits. In both panels, all SEDs have been normalized to unit flux at A lcst = 2500 A using linear interpolation between the adjacent data points (the 
fluxes are in units of erg s cirT 2 A -1 ). The J s —Ks > 2.3 selected objects have significantly redder rest-frame optical colours than the LBGs, hence we refer to 
them as "distant red galaxies" or DRGs. 
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FIG. 6. — Same as Figure|5]for the rest-frame SEDs of 2 < z < 3.5 DRGs in the FIRES fields. Only the Im galaxy template from Coleman et al. 1 1980) is 
shown here (black line) and used to illustrate the effects of interstellar extinction (grey lines) for Ay between and 3 mag, increasing in steps of 0. 3 mag. The Im 
template itself includes intrinsic reddening so that our simulations represent additional amounts of extinction. The extinction was applied using the Calz etti et al. I 
1 2000) law. 
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FIG. 7. — Diagnostic /gi47 s H K s colour diagrams to distinguish between dust obscuration and age effects, (a) H-K s versus / s —H. The diagonal line indicates 
the J s —K s = 2.3 mag colour cutoff used to select DRGs. The various track s show the evolution wit h redshift in the interval 2 < z < 4 of the observed colours of 
synthetic stellar populations. The template SEDs were generated with the Bruzual & Chariot i2003) synthesis code, assuming a Salpeter IMF between 0.1 and 
100 Mq and solar metallicity. The solid black lines correspond to models with constant star formation rate (CSF) and fixed ages of 10 Myr, 100 Myr, and 1 Gyr. 
The black dashed lines represent single stellar populations (SSPs) with fixed ages of 250 Myr and 1 Gyr. The black dotted line is the track for an SSP formed 
at zf = 5 and passively evolving with time. The grey solid lines show the effects of extinction on the 100 Myr old CSF model, using the Calzet tiet al. 1 1200(1) 
extinction law; the Ay increases from 1 to 6 mag (in steps of 1 mag) for the tracks from the bottom left to the top right, (b) Same as (a) for J s -H versus 7gi4— 7 S . 
(c) Colour distribution in the H—K s versus / s — H diagram of the combined sample of 2 < z, < 3.5 DRGs in the HDF-S and MS 1054-03 fields (black dots). The 
eiTor bars represent the 1<t measurements uncertainties and limits are plotted at the 2<r confidence levels. The shaded areas indicate the regions occupied by the 
tracks for dusty CSF models with fixed ages of 10 Myr (hatched region between dotted lines), 100 Myr (dark grey shaded region between solid lines), and 1 Gyr 
(light grey shaded region between dashed lines), (d) Same as (c) in the J s —H versus 7gi4 — J s diagram. In both diagrams, the spread of colours of DRGs in the 
direction perpendicular to the extinction paths is inconsistent with very young (< 100 Myr) populations highly obscured by dust; the presence of more evolved 
(> 250 Myr) populations with prominent Balmer/4000 A break is required for a large fraction of DRGs. 
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FIG. 8. — Colour excesses relative to dusty constant star formation models of the combined sample of 2 < z < 3.5 DRGs in the HDF-S and MS 1054 — 03 
fields, (a) Histogram of the deviation in J s — H colour of DRGs relative to the central locus of the dust-attenuated 100 Myr old CSF model tracks for 2 < z < 4 
in the H—K s versus 7 S — H diagram (based on FigureQ panels a and c). The one object with a limit on its J s — H colour is excluded. The black solid horizontal 
error bar indicates the maximum range in J s — H colour of the models around the central locus shown by the vertical black solid line at A(7 S — H) = 0. Similarly, 
the dotted and dashed grey error bars indicate the ranges for dusty CSF models of ages 10 Myr and 1 Gyr around their respective central locus, plotted with the 
corresponding vertical lines at their mean J s -H deviation from that of the 100 Myr old models, (b) Same as (a) for the deviation in Ig[4 — 7 S colour in the J s -H 
versus /gi4— 7 S diagram (based on Figure^ panels b and d). The reference set of models is again the suite of dusty 100 Myr old CSF populations, with central 
locus at A(/8|4 — 7 S ) = 0. The ten objects with limits on their 7gi4 -J s or J s — H colour are excluded. 
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FIG. 9. — Distributions of best-fit properties of the DRG samples in the FIRES fields. The panels to the left show the model results for a constant star formation 
rate (CSF) and to the right, for an exponentially declining star formation rate with e-folding timescale r = 300 Myr (T3ooMyr)- The models assumed a Salpeter 
IMF between 0. 1 and 1 00 Mq , solar metallicity, and the Calzetti et al. 1 2000 ) extinction law. The results are shown separately for the sample in HDF-S (hatched 
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FIG. 10. — Variations of the best-fit properties of DRGs with assumptions on the star formation history, extinction law, and metallicity. The median values 
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TABLE 1 

Summary of selected properties of the adopted / s - K s > 2.3 

SAMPLES AT 2 < Z < 3.5 IN THE HDF-S AND MS 1054-03 FIELDS 



HDF-S: 11 DRGs at Kf 1 < 22.5 mag 



Property Median Mean rms 



/Q ot (mag) 21.39 21.25 0.55 

V tot (mag) 26.33 26.69 0.96 

J s -K s (mag) 2.60 2.75 0.43 

z ... 2.50 2.60 0.42 



MS 1054-03: 23 DRGs at A r '° l < 21.7 mag 



Property Median Mean rms 



Xj 01 (mag) 20.95 20.84 0.74 

V tot (mag) 26.86 26.61 0.94 

J s -K s (mag) 2.79 2.85 0.43 

z ... 2.42 2.44 0.30 



NOTE. — All magnitudes are expressed in the Vega system. The 
■C and y tot data for the MS 1054-03 field sample are corrected for 
the (small) lensing magnifications. 



TABLE 2 

Summary of model results forJ s -K s > 2.3 samples in the HDF-S and MS 1054-03 fields" 



HDF-S 



Property 




Median 


CSF 
Mean 


rms 


Median 


T~300Myr 

Mean 


rms 


Age f sf 


(Gyr) 


1.7 


1.5 


0.86 


1.0 


1.1 


1.0 


SFR-weighted age (?}sfr 


(Gyr) 


0.85 


0.74 


0.43 


0.75 


0.93 


0.92 


Av 


(mag) 


2.7 


2.0 


0.78 


0.90 


1.20 


0.96 


M* 


(1O 1O M ) 


8.1 


15 


13 


7.5 


7.9 


5.7 


Lv.* b 


(10 10 Lv,s) 


5.0 


5.7 


3.6 


5.0 


5.7 


3.6 


Af*/£v,* 


(MqL^) 


1.2 


2.7 


2.2 


1.2 


1.5 


0.89 




(M yr-') 


91 


120 


96 


91 


150 


160 


SFR 


(M Q yr-') 


120 


150 


120 


6.7 


96 


140 


SFR/M* 


(10- 10 yr 1 ) 


7.6 


19 


26 


1.5 


18 


25 


SFR/Ly.* 


(10-'° Moyr-'L^) 


30 


33 


30 


1.2 


21 


32 



MS 1054-03 



Property 


Median 


CSF 
Mean 


rms 


Median 


T"300Myr 

Mean 


rms 


Age f sf 


(Gyr) 


2.0 


1.7 


0.84 


1.0 


1.3 


0.91 


SFR-weighted age (?}sfr 


(Gyr) 


1.0 


0.83 


0.42 


0.75 


1.1 


0.85 


Ay 


(mag) 


2.4 


2.2 


0.56 


0.90 


1.3 


0.80 


M* 


(10 10 M Q ) 


16 


21 


16 


10 


15 


11 




(10 10 Lv,s) 


5.1 


7.4 


6.4 


5.1 


7.4 


6.4 




(M L V ' ) 


2.3 


3.3 


1.9 


1.8 


2.4 


1.7 




(MQyr- 1 ) 


140 


150 


110 


94 


170 


150 


SFR 


(Movr" 1 ) 


170 


190 


130 


23 


69 


110 


SFR/M* 


(10- 10 yr 1 ) 


6.5 


15 


19 


1.5 


6.9 


10 


SFR/Lv.* 


(10- 10 Moyr-'L^g) 


31 


32 


22 


2.7 


11 


17 



" The samples selected with 2 < z < 3.5, with A" s < 22.5 mag for HDF-S and K s < 21.7 mag for the MS 1054 -03 field. The 
main m odel parameters are a Salpeter IMF between 0.1 and 100 Mq, solar metallicity, the extinction law of Calzetti et al. 
(2000), and either a star formation rate constant in time (CSF) or exponentially declining in time with e-folding timescale 
r = 300 Myr (r 3 ooM y r). 

b Th e Lv,± values were obtained independently, based on the less model-dependent method described by Rudnick et al. 
(2003), and are the dust-attenuated luminosities. These estimates are used in all ratios involving Lv,*. The ages, extinctions, 
stellar masses, and instantaneous SFRs were derived from the best-fit model to the SEDs. 



